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Abstract Phylogenetic relationships in the Dorcus velutinus species-group, con- 
sisting of D. velutinus, D. carinulatus, D. taiwanicus and D. japonicus, are examined based 
on the mitochondrial 16S rRNA gene sequences for re-evaluating the taxonomic status of 
Japanese D. japonicus and Taiwanese D. carinulatus. The relatively large sequence diver- 
gence in the D. velutinus species-group suggests that all the four taxa of this species-group 
(including D. japonicus and D. carinulatus) should deserve the status of valid species. The 
resultant phylogenetic trees indicate that the species-group is divided into two subgroups, 
one of which consists of D. velutinus and D. taiwanicus, whereas the other of D. japonicus 
and D. carinulatus. This agrees well with their hypothetical affinities based on their mor- 
phology. 


Introduction 


Dorcus velutinus J. THOMASON, 1862 was originally described from India, and 
later recorded also from Taiwan (Miwa, 1927) and Japan (TAKAI, 1981, misspelled as 
“D. veltinus" in his paper). NAKANE and MAKINO (1985) reexamined Taiwanese and 
Japanese specimens, and concluded that both the Taiwanese and Japanese forms were 
specifically distinct from D. velutinus from India. Consequently, they described the for- 
mer as Dorcus taiwanicus and the latter as D. japonicus (NAKANE & MAKINO, 1985). 
Besides these, other two species allied to D. velutinus have been described from Asia: 
i.e., D. ursulus ARROW, 1938 from Bhutan and D. carinulatus NAGEL, 1941 from Tai- 
wan. Of these, as mentioned below, D. carinulatus described on the basis of a single 
male specimen has been taxonomically problematical, because its diagnostic charac- 
ters mentioned in NAGEL’s (1941) description were not so clear, and, what is worse, the 
unique type specimen deposited in the Hannover Museum was destroyed during the 
World War II (KRAJCIK, 2003). BENESH (1960) suppressed D. carinulatus as a junior 
synonym of Metallactulus parvulus. MAES (1992) followed him, and when establishing 
a new genus, Velutinodorcus Maes, 1992 on the basis of D. velutinus as 1tS type 
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species, he moved three species allied to D. velutinus except D. carinulatus (i.e. D. ur- 
sulus, D. taiwanicus and D. japonicus) to this new genus. On the contrary, MIZUNUMA 
and NAGAI (1994) considered D. carinulatus to be a valid species, and regarded D. tai- 
wanicus as a junior synonym of the former. This arrangement was followed by ZHANG 
(1999) and KRAJICIK (2001, 2003). On the other hand, YANG and CHANG (1997) re- 
ported that D. japonicus (misspelled as *D. japonica" in their paper) was distributed in 
Taiwan as well as in Japan, and they illustrated its male genitalia. Afterwards, Li 
(2004) recognized D. japonicus and D. taiwanicus as two valid species occurring in 
Taiwan, but he made no comments on the taxonomic status of D. carinulatus. Quite re- 
cently, NAGAI and Fui (2005) revised D. velutinus and its related species including D. 
carinulatus under consideration. They concluded that D. carinulatus should be differ- 
ent from D. taiwanicus but identical with Taiwanese D. japonicus sensu YANG and 
CHANG (1997). Further, they regarded D. japonicus as a valid species. However, the 
taxonomic status of D. japonicus should be re-evaluated more carefully; according to 
YANG and CHANG 'S (1997) illustration, both the external and male-genitalic morphol- 
ogy of D. carinulatus (referred to Taiwanese D. japonicus by them) were quite similar 
to those of true Japanese D. japonicus. Thus, in the present study, we examine the phy- 
logenetic relationships among the members of the D. velutinus species-group (i.e., the 
members of the genus Velutinodorcus sensu MAEs’ (1992) and D. carinulatus) based 
on the mitochondrial DNA sequences to evaluate the taxonomic status of each species 
of this group including D. japonicus. 


Materials and Methods 


Insects 

We followed the taxonomic treatment of NAGAI and Fujii (2005). Four of the five 
species in the Dorcus velutinus species-group, D. velutinus, D. carinulatus, D. taiwani- 
cus and D. japonicus, were examined for DNA analysis. We selected Dorcus curvidens 
binodulosus and Dorcus titanus pilifer as outgroups for our analyses. See Table | for 
the detailed collection data of each beetle. 


Table 1. List of specimens used in the present study. 
eee 





Species Name Sampling Localities Accession No. 
Dorcus japonicus Isen, Tokuno Is., Kagoshima, Japan AB 236863 
Dorcus carinulatus Fengyuan, Taichung, Taiwan AB 236864 
Dorcus taiwanicus Habonshan, Nantou, Taiwan AB 236865 
Dorcus velutinus Doi Saket, Chiang Mai, Thailand AB 236866 
Dorcus curvidens binodulosus Yamanashi, Japan ABI 78 292* 
Dorcus titanus pilifer Sakyó, Kyoto, Japan AB178293* 


* Taken from previous studies (Hosoya et al.. 2001; HOSOYA & ARAYA, in press). 
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DNA extraction, amplification and sequencing 

We used specimens of adult or larva stocked in 99.5% ethanol. DNA was ex- 
tracted from small amounts of muscle using either the phenol/chloroform method 
(Hosova et al., 2001) or IsoQuickTM Nucleic Acid Extraction Kit (ORCA Research, 
Inc., USA). A part of the mitochondrial 16S rRNA gene was amplified by polymerase 
chain reaction (PCR) (SAiki et al., 1988) using six primers: 16SB. 16SC, 16SD, 16SH, 
16SK and 16SL (Hosoya et al., 2001; Hosoya & Araya, in press). The method of 
DNA amplification is described in detail elsewhere (Hosoya et al., 2001). Amplified 
DNA was purified by electrophoresis in 2% agarose gel, or PCR product pre-sequenc- 
ing kit (USB, Cleveland, USA). Nucleotide sequences were determined for both 
strands with a Thermo sequenase cycle sequencing kit (USB, Cleveland, USA) and a 
LI-COR Model 4200 Automated DNA Sequencer (LI-COR, Lincoln, USA), or 
SEQ2000 Dye Terminator Cycle Sequencing with Quick Start Kit and SEQ2000 
(Beckman Coulter, USA) using 16SB, 16SC, 16SD, 16SH and 16SH2 primers 
(Hosoya et al., 2001; HosovA & ARAYA, in press). 


Phylogenetic analysis 

Sequence data was aligned by CLUSTAL W 1.8 (THOMPSON ef al., 1994) at DNA 
Data Bank of Japan (DDBJ), using default gap penalties. The neighbor-joining (NJ) 
method was applied to infer the relationships among taxa based on pairwise distance 
matrix from KIMURA's two-parameter model (KIMURA, 1980), using CLUSTAL W. For 
maximum-parsimony (MP) analysis, branch-and-bound searches were performed 
using PAUP 4.0b10 (Sworrorb, 2002). The confidence levels of each branch were esti- 
mated using 1000 bootstrap replications (FELSENSTEIN, 1985). 

Baysian phylogenetic analyses were conducted with MrBayes 3.0 (HUELSENBECK 
& RONQUIST, 2001) under the general time reversible (GTR) model with proportion of 
invariable site and gamma shape parameter (GTR+I+G). Each Markov chain was 
started from a ramdom tree and run for 5X 10? generations, sampling the chains every 
100th cycle. All sample points prior to reaching stationary (1000 trees) were discarded 
as burn-in samples. Data remaining after discarding burn-in samples were used to gen- 
erate a majority rule consensus tree, where percentage of samples recovering any par- 
ticular clade represented the clade's posterior probability. 


Results 


Sequence variations 

All sequences are deposited in DDBJ nucleotide sequence databases under acces- 
sion numbers AB236863-AB236866 (Table 1). All sequences in the present study 
showed a substantial bias for adenine (A) and thymine (T) (average of 74.4%). Several 
observations have demonstrated the existence of a strong AT bias in the insect mito- 
chondrial genome (SIMON et al., 1994). 

The total mitochondrial 16S rRNA gene fragment consists of 996 bp, 176 (17.7%) 
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Table 2. Pairwise distance matrix from KIMURAS two-parameter model (KIMURA, 1980). 








Species Name | 2 3 + 5 
l. Dorcus japonicus 
2. Dorcus carinulatus 9.4 
3. Dorcus taiwanicus 10.2 10.1 
4. Dorcus velutinus 10.1 10.0 S3 
5. Dorcus curvidens binodulosus 16.6 18.4 16.6 16.6 
6. Dorcus titanus pilifer 23.8 23.4 21.3 21.4 20.7 





sites of which were variable. Aligned site of D. japonicus consists of 984 bp, D. cari- 
nulatus of 956 bp, D. taiwanicus of 977 bp, D. velutinus of 958 bp, D. curvidens bin- 
odulosus of 981 bp (data from Hosoya & Araya, in press, accession no. AB178292), 
and D. titanus pilifer of 976 bp (data from Hosoya et al., 2001, no. AB178293), re- 
spectively. 

The sequence divergence among taxa based on KIMURAS (1980) two-papameter 
model is summarized in Table 2. The sequence divergence within the D. velutinus 
species-group ranged from 5.3 (D. velutinus vs. D. taiwanicus) to 10.2% (D. taiwani- 
cus vs. D. japonicus), and those between the D. velutinus species-group and other con- 
geners using as outgroups varied from 16.6 (D. japonicus vs. D. curvidens 
binodulosus, D. taiwanicus vs. D. curvidens binodulosus, and D. velutinus vs. D. curvi- 
dens binodulosus) to 23.896 (D. japonicus vs. D. titanus pilifer). The sequence diver- 
gence between D. japonicus and D. carinulatus in problem are 9.4%, which is larger 
than those between D. velutinus and D. taiwanicus (5.3%). 


Phylogenetic relationships 

The NJ dendrogram derived from 16S rRNA sequences is shown in Fig. 1. The 
monophyly of the D. velutinus species-group was supported by high bootstrap propor- 
tion (BP) value (10096, node 1). Further, the species-group was divided into two clus- 
ters, both of which were supported by high BP value: one consisting of D. velutinus 
and D. taiwanicus supported by a BP value of 100% (node 2), and the other consisting 
of D. japonicus and D. carinulatus with a BP value of 9296 (node 3). 

As a result of MP analysis, a single most parsimonious tree was obtained 
(length=431, consistency index=0.85, retention index=0.54, rescaled consistency 
index=0.46). The resulting topology of the dendrogram of MP tree was identical with 
the NJ dendrogram, of which three nodes were supported by high BP valuez: 7096 
(100% for node 1, 78% for node 2 and 88% for node 3, respectively) (Fig. 1). 

The resulting topology of the dendrogram of Bayesian analysis was also identical 
with the NJ dendrogram, and posterior probabilities (PP) values showed that the three 
nodes of the Bayesian dendrogram were supported by high PP values20.90 (1.00 for 
node 1, 1.00 for node 2 and 0.93 for node 3, respectively) (Fig. 1). 
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Fig. 1. Neighbor-joining (NJ) dendrogram, with which the topology of those of MP and Bayesian analy- 
sis were identical, derived from Kimura’s (1980) two-parameter distances from the mitochondrial 16S 
rRNA sequences. Bar equals 0.05 substitutions/site. The numbers above internal branches represent 
the percentage of 1000 bootstrap for nodes supported (left, NJ; right, MP), and those below branches 
represent the Bayesian posterior probabilities. 


Discussion 


The sequence divergence in the mitochondrial 16S rRNA gene between D. japo- 
nicus and D. carinulatus (9.4%) observed in the present study is larger than that be- 
tween D. velutinus and D. taiwanicus (5.396), and as much as those between species in 
the other lucanid genus, Ceruchus (10.05%, between Japanese C. lignarius and Ameri- 
can C. piceus) (Hosoya et al., 2001). These data suggest that the difference in the se- 
quence divergence of mitochondrial 16S rRNA gene between D. japonicus and D. cari- 
nulatus should be considered as inter-specific level, and that the separation between 
them may have occurred earlier than that between D. velutinus and D. taiwanicus. 
These suggestions support the treatment of D. japonicus by NAGAI and Fuji (2005) for 
D. japonicus. 

As the results of present phylogenetic analyses, the monophyly of the D. velutinus 
species-group was supported by high bootstrap proportions (BP) value (100% in both 
of NJ and MP) and high posterior probabilities (PP) values (1.00). In addition, it was 
also revealed that the sequence divergence between each member of the D. velutinus 
species-group and the other Dorcus species used as outgroups was relatively large 
(16.6 to 23.4%). These results may suggest that Velutinodorcus Maes is possibly ap- 
propriate at least as one of the subgenera of the genus Dorcus, although further studies 
on the phylogenetic relationships among lucanid genera, including Gnaphaloryx 
BURMEISTER, 1847, as well as various species of the genus Dorcus are required. 
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The phylogenetic relationships within the D. velutinus species-group inferred in 
the present study indicate that this species-group is divided into two subgroups, one of 
which consists of D. velutinus and D. taiwanicus whereas the other of D. japonicus and 
D. carinulatus. This agrees well with their hypothetical affinities based on their mor- 
phological similarities (NAGAI & Fustl, 2005). The present results also revealed rela- 
tively remote phylogenetic relationship between D. taiwanicus and D. carinulatus, al- 
though both of them are distributed in Taiwan sympatrically (NAGAI & Fusu, 2005). 
Anyway, further detailed studies on the molecular phylogeny of the D. velutinus 
species-group including another species, D. urusulus, are strongly desired in order to 
discuss their divergence from the zoogeographical point of view. 
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z 約 


WEGE HII]: ミト コン ドリ ア 16S rRNA 遺伝 子 に 基づく Dorcus velutinus 種 群 の 分 子 
Rike, 日 4 よび 台湾 産 の 種 の 分 類 学 的 地位 に 関す る 検討 . 一 一 Dorcus carinulatus は , 
NAGEL (1941) に よっ て 台湾 中 部 か ら 得 られ た , た だ 1 個体 の 雄 に 基づい て 記載 され た が , タイ 
プ 標 本 が 第 二 verti 消失 し て し まっ た こと も あっ て , その 正体 が 不明 で , 記載 以降 , 538 
学 的 な 扱い が た び た び 変更 きれ て きた 問題 の 多い 種 で ある . 永井 ・ 藤 井 (2005)) は この の carinu- 
latus (NAGEL, 1941) を 含む 万 velutinus J. THOMASON, 1862 お よび その 近 緑 種 に 関し て 再 検討 し , 
楊 ・ 張 (1997) が 台湾 か ら か japonicus NAKANE et MAKINO, 1985 (ヤマ ト サ ビ クワ ガタ ) と し て 記 
$k L EZ ZD D. carinulatus (NAGEL, 1941) で あり , 台湾 に は D. taiwanicus NAKANE et MAKINO, 
1985 と D. carinulatus と の 2 種 が 分 布 し て いる と 結論 し た . この 見 解 は きま きわめ て 妥当 で ある が , 
D. carinulatus と D. japonicus は , 外部 形態 も また 雄 交尾 器 の : - ほ も よく 似 て お り , 両者 の 関係 に 

つい て は さら な る 精査 が 必要 と 判断 され た . そこ で 本 研究 で は , D. velutinus と その 近 緑 種 , D. 
carinulatus, D. taiwanicus, D. japonicus の 合計 4 種 に 関し て ミト コン ドリ ア 16S rRNA 遺伝 子 に 基 
づく 分 子 系 統 解 析 を 行い , 種 群 内 の 系 統 関係 を 明らか に する と と も に , D carinulatus & D. 
japonicus の 分 類 学 的 関係 に つい て も 考察 を 加え た . 解析 の 結果 , D. carinulatus © D. japonicus の 
間 の 塩基 置換 率 は 9.4% と , KEO D. velutinus と 台湾 産 の taiwanicus の 間 の 値 (5.4%) よ り も る 
高く , 既往 の ツヤ ハ ダ ク ワ ガ タ 属 の 日 本 産 (Cerzc/xs lignarius) と 北米 産 (C piceus) の 2 種 間 の 値 
(10.05%) に ほぼ 匹敵 する こと が 判明 し た . この こと か ら D carinulatus & D. japonicus は 遺伝 的 に 
も 別種 レベ ル の 違い が ある こと が 裏付け られ た . 解析 の 結果 得 ら れ た 系 統 関 係 か ら は , D. ve- 
lutinus TREE. D. velutinus © D. taiwanicus, B X OD. carinulatus と D. japonicus から な る 2 つの サ 
プ グ ルー プ に 分 れる こと が 判明 し た が , これ は 従来 の 外部 形態 に 基づく 本 種 群 内 の 類縁 関係 の 


Molecular Phylogeny of the Dorcus velutinus Group 529 


推定 結果 と も よく 一 致す る . また , D yes 種 群 全体 の 単 系 統 性 も きわ め て 高 し 確率 (NJ と 
MP 100% ; ベイ ズ の 事後 確率 1.00) で 支持 され た . Mass (1992) IX D. velutinus & を タイ プ 種 に Ve- 
lutinodorcus 属 を 創設 し た が , 今回 の 結果 は , この 見 解 に ある 程度 の 妥当 性 が ある こと を 示唆 し 
た も の と いえ る . 今後 は , 今回 の 解析 で は 使用 で き な か っ た の ursulus ARROW, 1938 は も ちろ ん , 
の ocs 属 の さま ざま な 種 や 他 属 の 種 も 加え た 系 統 解 析 を 行う こと で , D. velutinus 種 群 の 系 統 関 
係 を 明らか に する と と も に , Velutinodorcus 属 の 妥当 性 や D velutinus 種 群 の 系 統 生 物 地 理 に 関し 
て も 考察 を 加え て いき た い . 
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A Record of Pyrocoelia atripennis LEWIS (Coleoptera, Lampyridae) 
from Is. Okinawa-jima in the Ryukyus 


Masataka SATO 


DiaCuore Tokushige 306, Kamegahora 3-1404. Midori-ku, Nagoya, 458-0804 Japan 


Through the courtesy of Mr. M. Kimura and Miss A. OGarA, I was able to examine a spec- 
imen of Pyrocoelia atripennis LEWIS found in Is. Okinawa-jima. This was collected by an el- 
derly woman living in Okinawa. She said that the species sometimes came into the house during 
the winter. However, this species is naturally distributed in the Sakishima Islands. In this case, 
therefore, the firefly seems to have been brought in by someone from the Sakishimas, where it is 
very common in the winter. 

Specimen examined. 1 male, Kochinda, Is. Okinawa-jima, Ryukyus, 16-XI-2003、 H. 
KiNJO leg. 
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